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We present calculations of the Zeeman effect of narrow photoassociation lines of 40 Ca near the 3 Pi 
+ 1 So asymptote. Using a coupled-channel model we find a nonlinear Zeeman effect that even at 
low fields of a few mT amounts to several kHz. With this model we analyze previous measurements 
and give corrected long range dispersion coefficients of the 3 n u and 3 E+ states. 
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I. INTRODUCTION 

The scattering properties of alkaline earth atoms like 
Mg, Ca, Sr and the isoelectronic species Yb are cur¬ 
rently subject of intense investigations, both theoretically 
and experimentally. These properties are highly relevant 
e.g. for the generation of ultra-cold atoms in the quan¬ 
tum degenerate regime 0 - 0 , the evaluation of collisions 
that may limit the accuracy of the best optical clocks 
00, quantum computation 0, the determination of 
molecular potentials with high accuracy 0, or the pro¬ 
duction of ultra cold molecules 0 with the prospect of 
studying ultra-cold chemistry. Many of these investiga¬ 
tions are performed by using the photoassociation pro¬ 
cess, where a molecule is built from two ground state 
atoms in a resonant laser field. The intercombination 
transitions 1 So — 3 Pi of alkaline earth atoms with their 
extremely narrow natural line width of T/2-7r=7 kHz and 
T / 2 n = 374 ( 9 ) Hz for Sr and Ca, respectively, present al¬ 
most ideal model systems to perform such investigations 
with kilohertz accuracy [ff 11]. 

Despite of their similar electronic configuration, the 
twenty times smaller line width of Ca with respect to 
Sr leads to a striking difference for their long range po¬ 
tentials. In Ca the long range potentials are dominated 
by the van der Waals (Cq) interaction at the asymptote 
3 P + 3 S. Consequently, the small C3 coefficient (dipole- 
dipole interaction) in Ca results in two excited attractive 
potentials c 3 n u and a 3 E+ (Fig. [Qand Eq. |T]) which have 
a similar 1 /R e dependency as the ground state poten¬ 
tial [l 2 | ■ The c potential is purely attractive while the 
a potential with its Q = 1 component has a maximum 
of only h x 18 MHz at a distance of R = 5.5 nm due to 
the repulsive C3 term. The height of this maximum is 
small compared to the binding energy of the least bound 
state (« 983 MHz) which is still dominated by the Cq 
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FIG. 1: Relevant Hund’s case (a) molecular potentials of 
Ca 2 for photoassociation near the 3 Pi — 1 5'o asymptote (A = 
657 nm). The position of the transition is not to scale. 


coefficient. 

It was shown EH that the molecular states within these 
asymptotically similar potentials can couple by spin-orbit 
and rotational interaction remarkably even for the low¬ 
est angular momentum J = 1 . This behavior leads to 
a striking dependence of the molecular g factor on the 
rovibrational levels. The quadratic Zeeman effect that 
is significant in the photoassociation spectra of 88 Sr | 9 : ] 
is much smaller in Ca because of the different potential 
structure in Ca2 and Sr2 by the Cq and C3 terms, respec¬ 
tively. This results in a more widely spaced level struc¬ 
ture for Ca2 at the asymptote, whereas for Sr2 the levels 
of the fl = 0+ state (the plus sign referring to the total 
parity) are so close to the asymptote that the quadratic 
Zeeman effect is dominated by the coupling to the contin¬ 
uum above the asymptote. Thus the quadratic Zeeman 
contributions are significantly different in both molecules. 
In M it was therefore treated as being negligible under 
the experimental conditions and only considered in the 
uncertainty budget of the derived experimental g factors. 
In the present paper we give a more elaborate treatment 
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and calculate the Zeeman effect of the narrow photoasso¬ 
ciation lines of 40 Ca near the 3 Pi + 1 S 0 asymptote. Cor¬ 
recting previous measurements [lOJ we derive improved 
values of the Cq and C$ coefficients and g factors Hi¬ 
ll. ZEEMAN EFFECT OF AN ATOMIC PAIR 3 P 

+ x s 


The model is based on Hund’s case (a) potentials for 
the states a 3 £+ and c 3 II u (see Fig. [TJ. The short range 
potentials 0 are extended for R > 2 nm as 
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In 0 we calculated the molecular g factor from the 
expectation values of the Zeeman interaction of eigen¬ 
states of the field-free Hamiltonian; this corresponds to 
the first-order approximation. Here, we extend this the¬ 
oretical approach by diagonalization of the full Hamilto¬ 
nian being the sum of the field-free part and the Zeeman 
part. The first one is described in our previous work and 
the later one is set up as: 


Hz = HBgjBojz, ( 2 ) 


where gs is the Bohr magneton and gj the atomic g 
factor of state 3 Pj and j z is the projection of the total 
atomic angular momentum j onto the space fixed axis 
parallel to the external magnetic field with magnitude 
Bq. The actual Zeeman energy for magnetic fields be¬ 
low 10 mT is much smaller than the spin-orbit energy 
in Ca, thus recoupling of j to separate terms of orbital 
and spin angular momenta is unimportant. The atomic 
state 3 So does not contribute to the Zeeman energy. We 
neglect the Zeeman energy in the excited state by the 
pair rotation with the angular momentum l, for which 
the effective magnetic moment is only in the order of the 
nuclear magneton, and could give rise to shifts of not 
more than 1 kHz for 100 /iT field strength. 

The calculation is done in a Hund’s case (e) representa¬ 
tion | (L, S)j,l, J, M). The first three quantum numbers 
denote the atom 3 P ? , the constant quantum numbers for 
state 3 So are suppressed as well as ( LS ) in the state vec¬ 
tor below because of no importance here. J and M are 
the total angular momentum and its projection on the 
space fixed axis. The matrix elements of operator Hz 
(Eq. [2J are calculated with conventional angular mo¬ 
mentum algebra 


•w3j(./', 1, J, —M, 0, M) • w 6 j(j', J', l, J,j, 1) 
•[( 2 J+l)( 2 J' + l)( 2 j + l)(j + l)j ] 1 / 2 


( 3 ) 


w3j and w 6 j are Wigner’s 3j- and 6 j-symbols, respec¬ 
tively. The operator has non-vanishing matrix elements 




FIG. 2: Binding energies h ■ A b of the v’ = 1 lines in the 0 
and 1 potential. The red dashed line indicates the average of 
the M = +1 and M = — 1 energy. The dashed vertical line 
in the insets indicates the magnetic field of 285 /iT used in 
the experiment 0- In the insets the M — ±1 components 
are barely visible on the left side, because the linear Zeeman 
effect is huge compared to the quadratic one. 


for A J = 0, ±1 only. In the present experiments only low 
magnetic fields were applied, thus we truncate the ma¬ 
trix and span J < 3 for calculating the observed states, 
which have still fairly good quantum numbers J ~ 1. 

Fig. [2] shows the calculated Zeeman energy for the 
most weakly bound levels of state 0 + and l u demon¬ 
strating clearly the nonlinear behavior for magnetic fields 
below 10 mT. The excited ungerade state c0+ and the 
strongly mixed state (a,c)l u asymptotically approach the 
3 Pi+ 1 So state and correlate to the potentials denoted 
c 3 n u and a 3 £+ in Hund’s coupling case (a) shown in 
Fig. □ 

The inset in Fig. [2] gives a zoom at low magnetic field 
presenting the quadratic Zeeman shift of 30 kHz and - 
17 kHz for the M = 0 component of the least bound 
states in the two levels at the low field of 285 /iT used in 
the former experiment [10]. The dashed lines represent 
the average of the M = ±1 levels thereby canceling their 
linear Zeeman energy. At low magnetic field these aver¬ 
ages coincide well with the M = 0 level, which indicates 
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TABLE I: Measured and calculated binding energies h,A h (with J = 1) at zero field, their difference 5 and molecular g factors 
of the cOi and (a,c)lu states for vibrational level v' counted from the asymptote. Results of our calculation: Quadratic 
Zeeman effect described by scalar and tensorial susceptibility x scal and x tens , frequency shift Az], eman of the M — 0 component, 
calculated for a magnetic field of B = 0.285(7) mT. 


a small tensorial second-order Zeeman shift for the state 
l u . For comparison, the quadratic Zeeman shift of the 
atomic line l 1 Sq — 3 Pi), mostly resulting from admixture 
of the other atomic fine structure terms j = 0 and j = 2, 
is 6.3 x 10 7 Hz/T 2 [3] resulting in a shift of 5 Hz for the 
actual field and thus more than a factor of 10 3 smaller 
than for the weakly bound molecular states. 

The quadratic Zeeman shift was previously not in¬ 
cluded in the extrapolation of the measured energies to 
zero field 0- We calculated this shift A^ eemarl (col¬ 
umn 10 of Table QJ for the M = 0 component of all ob¬ 
served levels for 285 /j ,T and applied them as corrections 
to the data in 0 for obtaining the zero field frequencies 
Ap Xp (collected in column 4 of Table Q} . The experimen¬ 
tal molecular g-factors were derived from the M = ±1 
components correcting also the quadratic Zeeman effect 
are shown in column 7 of Table [I] Because the experi¬ 
ment did not distinguish between M = ±1 the sign of 
the g-factor is not determined. 

As in JTdJ the corresponding energies of the observed 
six PA resonances and their molecular g factors were fit¬ 
ted simultaneously. A non-linear least squares fit routine 
varies the slope of the short range repulsive branch (to 
get the proper phase of the long range wave function) and 
the long range parameters C$ and C$. The C^ coeffi¬ 
cient was held fixed at the value calculated from the life 
time of 3 Pi state of Ca, see ref. 0. The fit results for 
the energies and the g values are given in columns 3 and 6 
of Table Q] The corresponding long range parameters Cq 
and Cg are listed in Table HT1 together with theoretically 
derived values. Within the estimated uncertainty of 5% 
and 20% for the Cq and Cs parameters, respectively, our 
values agree with the theoretically determined ones. 

From the calculations we can extract the low-field coef¬ 
ficient of the quadratic Zeeman effect that can be written 
as a sum of a scalar and a tensorial part: 

A 4 2) (P, M)=-f [x scal + (3 M 2 - J(J + l)) X tens ] 

(4) 

Results are presented in Table [fl columns 8 and 9 fit- 


state 

Ce 

Cs 


10 7 cm 1 A 6 

10 9 cm -1 A 8 


this work 

HQ] 

\m 

[16] 

this work 

[10] 

m 

c 3 n u 

ffi 3 S u 

1.209 

1.335 

1.199 

1.353 

1.187 

1.313 

1.226 

1.358 

0.238 

0.905 

0.289 

0.818 

0.266 

1.057 


TABLE II: Long range parameters derived from the present 
work at the asymptote 3 Pj+ 1 So in comparison with theoret¬ 
ical results from Mitroy and Zhang [l6j] and Ciurylo et al. 

0 - 


ting the field range up to 0.3 mT of calculations as shown 
in Fig. [2] It is striking that the scalar part alternates 
in sign between the states Q = 0 + and fl = 1. This is 
directly related to the position of the perturbing states 
with respect to the observed states, which are all ,7=1 
with negative parity (noted below for short by a minus 
sign behind the J quantum number), i.e. levels with so- 
called e symmetry. From eq. [3] the non-diagonal matrix 
elements, which govern the quadratic Zeeman effect will 
have J' = 0 or J' = 2, thus the perturbing levels are 
/ levels with positive parity. We note that a quadratic 
contribution in B will also occur for A J = 0 because 
diagonal Zeeman terms will transform in non-diagonal 
terms (AQ = +/ — l)ina molecular basis because l and 
j from the atom pair basis are no longer good quantum 
numbers for the molecular case. Levels with J' = 0 and 2 
can only occur for states 0~, 1 and 2, where the first and 
last one correlate to the asymptotes 3 Po, 2 + 1 So and thus 
only accidental coincidences with the observed levels will 
give significant contributions. But the levels J' = 2— of 
Q = 1 correlating to the asymptote 3 Pi+ 1 So run parallel 
with the observed level structure and a regular pattern of 
the perturbation can be expected. We calculated these 
for the derived potential scheme and found that in all ob¬ 
served cases a level J' = 2—, Q = 1 is above the observed 
J = 1— level of Q = 0 + and below the one of Cl = 1. 
This leads to level shifts down in energy for the 0 + levels 
and up for 1 and with the sign definition of the scalar 
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susceptibility x scal in eq. [4] to positive and negative val¬ 
ues for 0+ and 1, respectively. The slow decrease of the 
magnitude of x scal for deeper bound levels corresponds 
to the increase of the level spacing going deeper into the 
potential. 

A similar calculation has been performed for Sr2 0. 
For the case 0 + the authors approximate the quadratic 
Zeeman shift as being inversely proportional to the dif¬ 
ference of the potential functions of 0 + and 1. This is 
justified because as mentioned in the introduction the 
quadratic Zeeman part for the bound levels very closely 
lying to the asymptote is dominated by the coupling to 
the continuum. Consequently, this approach fails for the 
observation of the deeper bound levels of state 1, because 
the perturbing level belongs to the same potential as the 
observed one and the energy difference originates from 
the rotational energy. This influence can be seen from 
significant deviations in Fig. 3 of reference [f|. 

III. CONCLUSION 

In conclusion we have calculated the influence of the 
quadratic Zeeman effect on the photoassociation lines of 


40 Ca near the 1 * 3 Pi + 4 5 6 7 8 So asymptote and shown that it 
is much smaller than in the case of Sr as a result of the 
very weak dipole - dipole interaction in calcium and the 
dominating van der Waals interaction in long range. The 
previous measurements have been corrected for this ef¬ 
fect since it was considerably larger than the estimated 
uncertainties for the unperturbed lines. Improved values 
of the long range parameters Cq and C$ have been de¬ 
rived that agree with the previously derived ones to a few 
percent and 20%, respectively, well within the previously 
estimated uncertainties. 
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